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The incipient infinite cluster appearing at the bond percolation threshold can be decomposed into 
singly- connected "links" and multiply-connected "blobs." Here we decompose blobs into objects 
known in graph theory as 3-blocks. A 3-block is a graph that cannot be separated into disconnected 
subgraphs by cutting the graph at 2 or fewer vertices. Clusters, blobs, and 3-blocks are special 
cases of /c-blocks with k = 1, 2, and 3, respectively. We study bond percolation clusters at the 
percolation threshold on 2-dimensional square lattices and 3-dimensional cubic lattices and, using 
Monte-Carlo simulations, determine the distribution of the sizes of the 3-blocks into which the blobs 
are decomposed. We find that the 3-blocks have fractal dimension d% — 1.2 ±0.1 in 2D and 1.1 5 ±0.1 
in 3D. These fractal dimensions are significantly smaller than the fractal dimensions of the blobs, 
making possible more efficient calculation of percolation properties. Additionally, the closeness of 
the estimated values for d% in 2D and 3D is consistent with the possibility that d% is dimension 
independent. Generalizing the concept of the backbone, we introduce the concept of a "/c-bone", 
which is the set of all points in a percolation system connected to k disjoint terminal points (or sets 
of disjoint terminal points) by k disjoint paths. We argue that the fractal dimension of a /c-bone is 
equal to the fractal dimension of /c-blocks, allowing us to discuss the relation between the fractal 
dimension of /c-blocks and recent work on path crossing probabilities. 



I. INTRODUCTION 

Percolation is the classic model for disordered systems 
]l|-f|. For concreteness we will study bond percolation 
systems in which bonds on a lattice are randomly occu- 
pied with probability p. Clusters are defined as groups of 
sites and bonds which are connected by occupied bonds. 
Clusters can be decomposed into objects known as blobs. 
Blobs are sets of sites and bonds which cannot be de- 
composed into disconnected sets by cutting only one 
bond. Equivalently blobs are sometimes described as be- 
ing multiply-connected — there are at least two disjoint 
paths between each point in a blob and every other point 
in the blob. The decomposition of the entire percolation 
cluster into blobs has been extensively studied Q , as has 
the distribution of sizes of blobs in the backbone [|| . For 
both cluster and backbone blobs, the fractal dimension 
of the blobs is the fractal dimension of the backbone. 

Here we address the questions of (i) whether there are 
more fundamental objects into which blobs can be de- 
composed, and (ii) whether these objects then be further 
decomposed. To answer these questions, we employ the 
language of graph theory, in which sites are the vertices 
and bonds are the edges of a graph ||. 

One can define /c-connected graphs (or /c-blocks) as 
graphs which cannot be separated into disconnected sub- 
graphs by cutting the graph at fewer than k vertices . 
Thus, clusters are 1-blocks and blobs are 2-blocks. The 



natural next level of decomposition of percolation sys- 
tems is to decompose blobs(2-blocks) into 3-blocks. By 
the definition above, 3-blocks are graphs which cannot 
be decomposed by cutting the graphs at fewer than 3 
vertices. From a physicist's point of view, one can un- 
derstand what 3-blocks are by considering a blob as a 
resistor network with each bond being a resistor. As- 
sume one is trying to determine the resistance between 
2 vertices of the network. One can simplify the network 
by using Kirchoff's Laws to replace groups of sequen- 
tial bonds and groups of parallel bonds by single virtual 
bonds having resistance equivalent to the bonds replaced. 
After this has been done as completely as possible, what 
are left are 3-blocks. We define the mass of a 3-block 
as the number of virtual bonds plus the number of non- 
replaced original bonds remaining in the 3-block. Figures 
|l| and H provide examples of the decomposition of a blob 
into 3-blocks. It has been shown || that the decomposi- 
tion of 2-blocks into 3-blocks is unique. 

Determining the scaling of the distribution of the 3- 
blocks into which the 2-blocks can be decomposed is the 
subject of this paper. In graph theory, the sites are typ- 
ically not constrained to a lattice structure, and one is 
only concerned with the topology of the graphs; we will, 
however, work on square and cubic lattices. 
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II. NOTATION 

Because we deal with a number of different types of 
fractal objects, we employ the following notation: 

(i) The fractal dimension of an object of type X will 
be denoted as dx- 

(ii) The number distribution of objects of type X in 
space Y of size L will be denoted as n(7Vx, Ly). 

(iii) The exponent of the power-law regime of a distri- 
bution of objects of type X in space of type Y will 
be denoted as tx,y- 

(iv) The amplitude of a distribution of objects of type 
X in space of type Y will be denoted as Ax,y- 

(v) We define d n Y through the relation 

(n(L)) ~ L dnY , (1) 

where (n(L)) is the average number of disjoint ob- 
jects of a given type in space Y. 

(vi) We use 0,1,2,3... to denote /c-blocks with k = 
0, 1, 2, 3 . . . corresponding to Euclidean space, clus- 
ters, blobs, and 3-blocks respectively. We use B to 
denote the percolation cluster backbone. 

(vii) Additionally, because, as noted above, objects such 
as 3-blocks can be nested, we denote quantities 
that relate to all levels of nesting with an asterisk. 
Specifically, r x Y and A x Y denote the exponent of 
the power-law regime and the amplitude of a dis- 
tribution of nested objects of type X at all levels of 
nesting in space of type Y. Similarly, d* y is defined 
through the relation 

{n*(L))~L d nr, (2) 

where (n*(L)) is the average number of nested ob- 
jects at all levels of nesting of a given type in space 
Y. Quantities not qualified with an asterisk will de- 
note quantities at a single level or quantities which 
cannot be nested. 

Using this notation, previous results are || 

n(N 2 , L B ) = A 2 , B L d ~* N~ T2 - B f L (j£) (3) 

for the number distribution of blobs of mass N2 in the 
percolation cluster backbone and || 

n(JV 2 ,L!) = A^L^N-^h (j£) (4) 

for the number distribution of blobs of mass N2 in the 
whole percolation cluster. The finite-size scaling function 



m Eqs. (Q) and ([|) approaches when x > 1 and 
is 1 otherwise. 

In analogy with Eqs. (||) and @) we expect the number 
distribution of 3-blocks at all levels of nesting in a blob 
to be 

n*(N 3 ,L 2 ) =Al 2 L<*N- T ^f c (^-) f L (j£) , (5) 

where c is the mass of the smallest 3-block and the finite- 
size scaling function f c (x) approaches when x < 1 and 
is 1 otherwise, reflecting the fact that there cannot be 
any 3-blocks smaller than the smallest size c. In all di- 
mensions and for all lattices, c = 5. For simplicity we 
will approximate n*(AT 3 , Lb) as 

n*{N 3 ,L 2 ) = {Al 2 L<*N- Tl > c < N 3 < aL d * (6) 
1 otherwise. 



III. SIMULATIONS 

We perform simulations with p = 0.5, the exact perco- 
lation threshold for 2D [g,g andp = 0.2488126, the most 
precise current estimate for the percolation threshold for 
3D ||. We created percolation clusters which included 
the sites (0, L/2) and (L, L/2) for the 2D simulations and 
the sites (0, L/2, L/2) and (L, L/2, L/2) for the 3D simu- 
lations, decomposed the backbones determined by these 
sites into blobs and then decomposed the blobs into 3- 
blocks. We study both distributions of 3-blocks in blobs 
of given mass, iV~2, and distributions of 3-blocks in back- 
bones in systems of a given size, L. For purposes of analy- 
sis, we group together blobs with mass 2 m_1 < N2 < 2 m . 

We perform the decomposition into 3-blocks along the 
lines of the procedure sketched in Ref. ||. Basically, 
this procedure is as follows: We first designate the blob 
that we are decomposing as the 2-block graph G. The 
natural next level of decomposition is to identify con- 
nected subgraphs with two or more edges that are con- 
nected to G at only two vertices. We denote these sub- 
graphs Gi, G2, G3, . . . of G as 2-terminal objects. These 
2-terminal objects can then be replaced in G by "virtual 
edges," ei, e2, e3, . . . Note that this process can be con- 
tinued recursively. That is, the subgraph Gi may itself 
contain sub-graphs, G^i, G$2, G^3, . . . that are connected 
to Gi at only two vertices; we then replace the subgraphs 
Gij in Gi by virtual edges eoij- The process continues 
until the only remaining subgraphs are those that cannot 
be decomposed further by making cuts at two vertices; 
these, by definition, are 3-blocks. An example of this 
decomposition is shown in Fig. ^ Other methods of de- 
compostion into 3-blocks are described in Refs. [p|]io|]. 

We perform at least 3700 realizations for each system 
size; for the smaller system sizes for which the simula- 
tions run more quickly we performed as many as 10 8 re- 
alizations. Because, the larger the systems the larger the 
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number of 3-blocks contained in the system, the statistics 
for the larger systems was acceptable despite the lower 
number of realizations. We bin the results for all system 
sizes in order to smooth the plots. 



IV. TWO SPATIAL DIMENSIONS 

In this section we discuss our results for 3-blocks in 
2D percolation. Results in 3D are analogous and are dis- 
cussed in the next section. 



A. 3-blocks in Blobs 

Figure |(a) plots the distributions P*(N 3 \N 2 ), the 
probability that a 3-block contained in a blob of size N 2 
contains Ns bonds, for various values of N 2 . P(Ns\N 2 ) is 
the number distribution n*(A^3, A^) normalized to unity. 
Consistent with Eqs. (||) and @, the plots exhibit power- 
law regimes followed by cut-offs due to the finite size of 
the blobs. The "bumps" in the distributions right be- 
fore the cutoffs represent 3-blocks which would have been 
larger but are truncated due to the finite size of the blobs 
in which they are embedded. We estimate the slope of 
the power law regimes, 73 2 , to be 2.35 ± 0.05. Since 



and 



we expect 



7V"o no L 



N 2 ~ L d 



d 3 



Ns 



N d 3 /d 2 



(7) 
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In Fig. ||(b), we show the collapsed plots in which we 
scale the distributions by N 2 ds ^ d2 using the most precise 
published estimate for d 2 , 1.6432 ± 0.0008 Q. (A con- 
sistent more recent estimate, d 2 = 1.6431 ± 0.0006, is 
given in Ref. [Q.) Visually, we find the best collapse is 
obtained for d 3 = 1.20 ± 0.1. 



We can also estimate ds using Eq. dA5|) from the ap- 
pendix 



3,2 



1) 



a n2 



d 2 



(10) 



Using r 3 * 2 = 2.35 ± 0.05 and d 2 = 1.6432 ± 0.0008, 
results in an estimate of ds = 1.22 ± 0.05. 



B. 3-blocks in Backbone 

Figure ||(a) plots the distributions P*(7V 3 |L jB ), the 
probability that a 3-block contained in the backbone of 
a system of size L contains A/3 bonds, for various values 
of L. P*(N 3 \L B ) is the number distribution n*(Ns,L B ) 



normalized to unity. Consistent with Eqs. (||) and (^), 
the plots exhibit power-law regimes followed by cut-offs 
due to the finite size of the systems. We estimate the 
slope of the power law regimes, t£ b , to be 2.25 ± 0.05. 
In Fig. ^(b), we show the collapsed plots in which we 
scale the distributions by L ds . Visually, we find the best 
collapse is obtained for ds = 1.15 ± 0.1. 

Next we consider the distribution of "top-level" 3- 
blocks in the backbone. Top-level 3-blocks are those not 
contained within another 3-block. In Fig. ||(a), we plot 
the distributions P(Ns\Lb), the probability that a top- 
level 3-block contained in the backbone of a system of 
size L contains Ns bonds, for various values of L. The 
plots exhibit power-law regimes followed by cut-offs due 
to the finite size of the systems. The exponent of the 
power-law regimes t 3 ,b is estimated to be 1.6 ± 0.05. In 
Fig. ||(b), we show the collapsed plots, in which we scale 
the distributions by L ds . The best collapse is obtained 
for ds = 1.15 ±0.1, the same value as for the distributions 
of 3-blocks of all levels. Thus the fractal dimensions of 
the top level 3-blocks is the same as the fractal dimension 
of 3-blocks of all levels but the slopes of the jpower law 
regimes are different; th is is seen also in Fig. 

We can also use Eq. ( |A10| ) 



ds(rs,B - 1) = d nB -. 

v 



(11) 



to obtain an estimate of ds. Since d n B is known ex- 
actly in 2 dimensions and has been well studied in higher 
dimensions and because one can usually determine the 
slope ts,b more accurately than ds can be determined by 
finding the best scaling collapse, we determine ds more 
accurately by solving Eq. (|l|) for ds- Using our esti- 
mate for ts,b above we find ds = 1.25 ± 0.1. Combining 
this result with our earlier estimates, we make the final 
estimate 



ds 



1.20 ±0.1. 



(12) 



C. Why the Fractal Dimension of 3-blocks is Smaller 
than the Fractal Dimension of the Backbone and 
2-blocks 

The fractal dimension of the 3-blocks is considerably 
smaller than the fractal dimension, ds = 1.6432 ±0.0008 
iPH , of 2-blocks (blobs). This is because virtual bonds 
are counted as one bond even though they replace many 
bonds. This can be seen if we plot the distributions 
P*(Ms\Lb), the probability that a 3-block contained in 
the backbone of a system of size L contains Ms bonds 
where we can count not the virtual bonds, but all bonds 
contained in a 3-block. In Fig. 0(a) we plot P*(M 3 \Lb) 
for various L. The best collapse for these plots (Fig. 0(b)) 
corresponds to a fractal dimension of 1.6 ± 0.1 consistent 
with the fractal dimension of 2-blocks in 2D. This can be 
understood as a reflection of the fact that in a system of 
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size L, the mass of the largest 3-block(counting all bonds) 
can be the same as the backbone mass. This is similar to 
the situation with blobs and backbones; the largest blob 
in a backbone can be as large as the whole backbone, 
which explains why the fractal dimension of blobs is the 
same as the fractal dimension of the backbone. 

Replacing a group of bonds by a virtual bond is anal- 
ogous to removing dangling ends on a cluster when de- 
termining the backbone. 

V. THREE SPATIAL DIMENSIONS 

Our analysis of the results of the 3D simulations pro- 
ceeds in a similar manner to the analysis for 2D. 

A. 3-blocks in Blobs 

Figure |(a) plots the distributions P*(iV 3 |7V 2 ), the 
probability that a 3-block contained in a blob of size N 2 
contains Ns bonds, for various values of N 2 . We estimate 
the slope of the power law regimes, r| 2 , to be 2.63 ±0.05. 
In Fig. ||(b), we show the collapsed plots in which we 
scale the distributions by N 2 ds/d2 with d 2 = 1.87 ± 0.03 
. Visually, we find the best collapse is obtained for 
d 3 = 1.15 ±0.1. 

Estimating ds using Eq. ( |A5[ ) from the appendix 

d3(>3* 2 - 1) = < 2 = d 2- (13) 

with r 3 * 2 = 2.63 ± 0.05 and d 2 = 1.87 ± 0.03, results in 
an estimate of ds — 1.15 ± 0.05. 



B. 3-blocks in Backbone 

Figure ^|(a) plots the distributions P*{N%\Lb), the 
probability that a 3-block contained in the backbone of a 
system of size L contains A/3 bonds, for various values of 
L. We estimate the slope of the power law regimes, T3 B , 
to be 2.55±0.05. In Fig. ^(b), we show the collapsed plots 
in which we scale the distributions by L ds . Visually, we 
find the best collapse is obtained for ds = 1.15 ± 0.1. 

Next we consider the distribution of "top-level" 3- 
blocks in the backbone. In Fig. |l0|(a), we plot the dis- 
tributions P(Ns\Lb)^ the probability that a top-level 3- 
block contained in the backbone of a system of size L 
contains Ns bonds, for various values of L. The expo- 
nent of the power-law regimes ts,b is estimated to be 
2.0 ± 0.05. In Fig. |l0|(b), we show the collapsed plots, in 
which we scale the distributions by L ds . The best col- 
lapse is obtained for ds = 1.15 ±0.1, the same value as for 
the distributions of 3-blocks of all levels. As in 2D, the 
fractal dimensions of the top level 3-blocks is the same 
as the fractal dimension of 3-blocks of all levels but the 
slopes of the power law regimes are different; this is seen 
also in Fig. |ll|. 



Using Eq. ( |aTc| ) 

d 3 (T 3 ,B ~ I) = d nB = J. (14) 

to obtain an estimate of ds with our estimate for ts,b 
above we find ds = 1.14 ± 0.1. Combining this result 
with our earlier estimates, we make the final estimate 

d 3 = 1.15 ±0.1. (15) 

The simulation results notwithstanding, it would be 
surprising if ds were smaller in 3D than in 2D because, 
below the critical dimension d c = 6, both the fractal 
dimensions of clusters and blobs increase with the Eu- 
clidean dimension. This suggests that while the actual 
values of ds may be within the bounds we have estimated, 
the actual values will be consistent with ds (2D) < 
ds (3D). 

As in 2D, if we do not replace two-terminal objects in 
a 3-block by a single virtual bond, the fractal dimension 
of the 3-block is that of a blob(see Fig. fL2|) . 

Estimates for all of the 2D and 3D exponents are sum- 
marized in Table I. 



VI. DECOMPOSITION OF THE WHOLE 
PERCOLATION CLUSTER 

While we have only decomposed 2-blocks that com- 
prise the cluster backbone, we could proceed similarly 
for all 2-blocks into which a cluster is decomposed. The 
fractal dimension of the 3-blocks into which a cluster is 
ultimately decomposed should be the same as the fractal 
dimension of the 3-blocks into which the backbone is ulti- 
mately decomposed. The only difference we would expect 
in our results would be that the slope of the power-law 
regime of the distribution of top-level 3-blocks would be 
given by 

^3(^3,1 - 1) = d n i = d. (16) 
the analogy of Eq. (|A10[) 

VII. K-BONES AND PATH CROSSING 
PROBABILITIES 

Just as blobs and backbones have the same fractal di- 
mension, we can identify objects analagous to backbones 
which have the same fractal dimensions as /c-blocks. We 
define a /c-bone as the set of all points in a percolation 
system connected to k disjoint terminal points(or sets of 
disjoint terminal points) by k disjoint paths. Thus the 
backbone is a /c-bone with k = 2. Just as the largest 
fc-blocks into which a backbone can be decomposed are 
2-blocks, the largest /c-blocks into which a /c-bone can 
be decomposed are /c-blocks. The fractal dimension of 
/c-bones is the fractal dimension of the /c-blocks. One 
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can see this easily by noting that if the k terminal points 
which define a /c-bone are connected to each other, the 
resulting structure is /c-block. 

Recent work |Q has identified a relationship between 
path crossing probabilities and the fractal dimensions of 
percolation structures. Specifically, consider the proba- 
bility, that in an annular region the small inner circle 
of radius r is connected to the larger outer circle of radius 
R ^> r, by k disjoint paths. Then 

an 

It has been observed Q that x \ is the codimension of 
the percolation cluster and £2 is the codimesion of the 
backbone. We extend these observations to the case of 
general k 

d-x k = d k , (18) 

where d is the spatial dimension of the system. This 
should hold in all dimensions where the annulus is now 
defined by two hyper spheres. It has been argued [|l2| that 

X k < X k < x 2k , (19) 

where x k is the polychromatic path crossing exponent 
|L4| and which has been found rigorously in 2D to be iQ 

Xk = T2 (fc2 " 1} - (20) 

Using Eqs. (|l|), @, and @, we find in 2D 

- s < * < 5- < 21 > 

consistent with our estimate for ds in 2D. 

The relationship between the path crossing problem for 
k = 2 and the backbone dimension has been recently ex- 
ploited to determine ds very accurately using a transfer 
matrix technique Possibly similar methods can be 

employed to find the fractal dimension of fc-bones (and 
therefore ^-blocks) with k > 3 to high precision. 

VIII. RELATIONSHIP TO RENORMALIZATION 
GROUP 

The process of replacing a 2-terminal object, t, by a 
single virtual bond and then replacing 2-terminal objects 
within t by single virtual bonds and so on is reminiscent 
of the decimation process in renormalization group (RG) 
approaches to percolation ||,||,|l5],[L6| . It is here, how- 
ever, that the similarity ends. The decimation process 
performed in the decomposition into 3-blocks is an exact 
decimation performed on objects in individual realiza- 
tions while the RG decimation is performed on the lattice 
and is an approximation, except for hierarchical lattices. 
Also, the purpose of the decomposition into 3-blocks is 



to improve computational performance and analyze the 
properties of substructures of the cluster while the pur- 
pose of RG calculations is to find properties of perco- 
lation analytically. Finally, whereas RG approaches on 
hierarchical lattices result in objects that are finitely ram- 
ified, the decomposition into 3-blocks we perform main- 
tains the infinite ramification of the Euclidean lattice. 



IX. COMPUTATIONAL IMPLICATIONS 



The fact that the fractal dimension of 3-blocks is sig- 
nificantly smaller than the fractal dimension of 2-blocks 
has important computational implications. We can effi- 
ciently calculate properties(e.g. resistance, velocity dis- 
tributions, self avoiding walk statistics) of a percolation 
cluster or backbone as follows: 

(i) decompose the cluster or backbone into 2-blocks 

(ii) decompose the 2-blocks into 3-blocks 

(iii) calculate the desired properties of the 3-blocks 

(iv) algebraically determine the properties of the 2- 
blocks from the properties of the 3-blocks 

(v) algebraically determine the properties of the cluster 
or backbone from the properties of the 2-blocks. 

In many cases the computation will require less CPU 
(computer processing) resource when the complexity of 
the computation is a power law or exponential of the 
mass of the object for which the property is being cal- 
culated. By decomposition we make the mass of these 
objects smaller. Reduced CPU resource usage is also 
obtained if only a decomposition into 2-blocks is made 
although the saving are less. Systems of larger size than 
could be treated before can now be treated when we de- 
compose into 3-blocks because the fractal dimension of 
the 3-blocks are lower than that of the object in which 
they are embedded; this is not true if we only decompose 
into 2-blocks. 

As an example of the dramatically smaller size of the 
largest 3-block versus the size of the largest blob con- 
sider a 3D system of size L = 1000. At criticality, the 
largest mass blob in the backbone will be of the order 
£i.62 ^ 63,000 while the mass of the largest 3-block in 
the backbone will be of L 1 ' 2 ~ 4000. In Fig. || we show 
an actual simulation realization in which a blob of 950 
bonds is decomposed into a 3-block with only 216 virtual 
bonds, greatly reducing the computational complexity. 
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X. DISCUSSION 

Traditionally the decomposition of percolation systems 
has been to decompose the system into clusters (1-blocks) 
and to decompose the clusters into blobs (2-blocks). We 
extend this decomposition by decomposing 2-blocks into 
3-blocks. 3-blocks are especially interesting because in 
contrast to 1- and 2-blocks, the 3-blocks have the prop- 
erty that they can be nested. That is, 2-terminal objects, 
which are replaced by single virtual bonds in a 3-block, 
can themselves contain other 3-blocks. Because of this 
replacement of a 2-terminal object by a virtual bond, 
the fractal dimension of 3-blocks is significantly smaller 
than the fractal dimension of 2-blocks. As discussed in 
the previous section, this smaller fractal dimension has 
important computational implications for the size of per- 
colation systems which can be analyzed and the speed at 
which the analysis can be performed. 

In addition, within the error bars of our calculations, 
the values for the 3-block fractal dimension appear to be 
identical for 2D and 3D systems. Simulations of larger 
systems and higher dimension systems could help answer 
whether in fact d% is independent of dimension ("super- 
universal"). It will also be of interest to determine the 
properties of /c-blocks with k > 3. 
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APPENDIX A: RELATIONSHIPS AMONG 
EXPONENTS 

Here we ask if any of the fractal dimensions and power- 
law regime exponents we have identified are related. To 
answer this question we must first briefly review some 
existing results for relations between other exponents. 



1. Previous Results 

It has been shown generally 0,(l8| that, for disjoint 
objects of type X embedded in a space Y, 



dx(rx,Y - 1) = d nY . 



(Al) 



Equation (|Al[) holds if tx < 2 or if d n y is equal to the 
fractal dimension of space Y, dy. 

Special cases of Eq. ( |Al| ) have been identified previ- 
ously for Y — 0, 1, and 2 corresponding to Euclidean 
space, percolation cluster space and percolation back- 
bone space, respectively. 



(i) The first is the familiar scaling relation for the 
Fisher exponent r ||[| 



df(r — 1) = d 



n0 



d 



clusters, (A2) 



where d is the Euclidean dimension, df the frac- 
tal dimension of the cluster, and r the exponent of 
the power-law regime in the distribution of cluster 
sizes. 

(ii) In Ref. Q it was shown that 

^biob-ci(rbiob-ci - 1) = d n i = df, (cluster blobs) 

(A3) 

where (ibiob-ci and Tbiob-ci are the fractal dimension 
and the exponent of the power-law regime, respec- 
tively, for all blobs in the cluster. 

(iii) In Ref. Q it was argued that 

dbiob-bb(7biob-bb - 1) = d nB = dred, (backbone blobs) 

(A4) 

where dbiob-bb and Tbiob-bb are the fractal dimension 
and the exponent of the power-law regime, respec- 
tively, for those blobs in the backbone and d re( ± is 
the fractal dimension of singly-connected red bonds 
in the backbone. 

Both (ibiob-ci and <ibiob-bb are equal to ds, the back- 
bone fractal dimension. In (i) and (ii), Eq. (|A1[) applies 
because d n Y = dy] in (iii), Eq. (Al) applies because 
rx < 2. 



2. 3-blocks in Blobs 

In analogy with Eqs. ( |A2| ) and (|A3|), we would expect 



*(r 3 * 2 -l) =< 2 = d 2 . 



(A5) 



We first confirm that the total number of 3-blocks in 
blobs scales with the exponent d 2 . If 



< n(L) >r 



and 



N 2 ~ L d 



(A6) 



(A7) 



then we would expect 



< n(N 2 ) >~L dn2/d \ (A8) 

Figures |l^(a) and [l3](b) are log- log plots of (n(N 2 )), the 
average number of all 3-blocks in a blob, versus blob size 
N 2 for 2D and 3D, respectively. The straight line fits 
with slope 1.0 ± 0.05 are consistent with d n2 = d 2 . Our 
simulation results in 2D from Section IV, ds = 1.20 and 
?3,2 = 2.35 result in ds(r^ 2 — 1) = 1.62 close to the value 
d 2 = 1.6432. In 3D, our simulation results from Section 
V, d 3 = 1.15 and r 3 , 2 = 2.63 result in d 3 (r 3 * 2 - 1) = 1.87 
identical to the value d 2 = 1.87 Ol. 
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3. 3-blocks in Backbone 

Because the number of top-level 3-blocks in the back- 
bone is proportional to the number of 2-blocks in the 
backbone, the number of top level 3- blocks in the back- 
bone should scale the same way the number of 2-blocks 
in the backbone. For all dimensions and lattices, d n B has 
been shown to be |l9| , ^Q[| 

d n B = d red ^, (A9) 

where v is the exponent associated with the divergence 
of the correlation l ength as p approaches p c [1,2]. In 2D 
\jv is exactly 3/4 felLb2[ ; in 3D, \ jv has been estimated 
to be 1.143 ± 0.01 f^lj. We would expect 

d 3 (T 3 ,B -1) =d nB = -. (A10) 

Figures |l4|(a) and |l4|(b) are log- log plots of (n(L)), the 
average number of top-level 3-blocks in the backbone, 
versus system size L for 2D and 3D, respectively. The 
straight line fits with slope 0.75 ±0.05 and 1.14 ±0.05 are 
consistent with the exact and previously estimated val- 
ues for \ jv of 3/4 and 1.143 in 2D and 3D, respectively. 
Our 2D simulation results from Section IV, ds = 1.15 and 
73,2 = 1.60 result in ds(rs,B — 1) = 0.69 close to the value 
\jv = 3/4. For 3D, our simulation results from Section 
V, ds = 1.15 and rs^ — 2.0 result in ds(rs y B — 1) = 1.15 
close to the value 1/v = 1.143. 
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TABLE I. Measured fractal dimension, measured power-law regime exponent, and calculated fractal dimension for 3-blocks 
in 2D and 3D. The calculated value of ds is determined by Eq. ( |lo| ) for 3-blocks in a blob and Eq. ( [ill ) f° r 3-blocks in the 
backbone. 



2D 












T 






ME A STIR ED 

1V± J_j Jr\. O KJ A. VJ_jJ_-/ 


MF, A STIR ED 

1V± J_j Jr\. O KJ ±VJ_jJ_-/ 


n ALCUT ATED 


All ^ KlonlrG in KloK 
.ri.ll O Ul<J^ls.o 111 U1UU 


1.20 ± 0.1 


O O K _L 
^j.kJO _l_ \J.\JO 




All 3 blocks in backbone 


1.15 ±0.1 


2.25 ±0.05 




Top level 3 blocks in backbone 


1.15 ±0.1 


1.60 ±0.05 


1.25 ±0.1 


3D 










d 3 


T 






MEASURED 


MEASURED 


CALCULATED 


All 3 blocks in blob 


1.15 ±0.1 


2.63 ±0.05 
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1.15 ±0.1 


2.55 ±0.05 




Top level 3 blocks in backbone 


1.15 ±0.1 


2.0 ±0.05 


1.14 ±0.1 
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FIG. 1. (a) Decomposition of 2-block G into subgraphs Gi, Gi, and C3. The rightmost graph represents G with the 
sub-graphs replaced by equivalent "virtual edges." (b) Subgraph G2 of G is decomposed by identifying subgraph G21. The 
rightmost graph represents G2 with the subgraph G21 replaced by it equivalent edge, (c) Subgraph G3 of G is decomposed by 
identifying subgraph G31. The rightmost graph represents G3 with the subgraph G31 replaced by it equivalent edge. In (a), 
(b), and (c) virtual edges are denoted by dashed lines. Note that while not shown in this figure, subgraph G31 could be further 
decomposed. The 3-blocks contained in the graph G are G21, having 5 edges, and G3 (with the subgraph G31 replaced by its 
equivalent edge) having 8 edges. 
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FIG. 2. Example of decomposition of backbone into 3-blocks. The thin lines represent the bonds in the backbone between 
points {0,15} and {31,15} on a lattice with L=32. The backbone is composed of a few single bond blobs connected to the 
terminal points and a single large blob containing 950 bonds. The thick lines represent the virtual bonds of a single top-level 
3-block into which the blob has been decomposed. This 3-block contains 216 virtual bonds. Some of the groups of bonds 
replaced by virtual bonds can themselves be decomposed into lower level 3-blocks and so on. 
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x=N 3 /N 2 d3/d2 

FIG. 3. 2D (a) Distributions P*(N3\N2) of the number of 3-blocks of mass N3 in a blob of size N2 versus N3 for (from bottom 
to top) N2 = 2 10 , 2 12 , 2 14 , and 2 16 . The distributions exhibit a power- law regime with slope —2.35 =b 0.05 (b) Distributions for 
N2 = 2 12 , 2 13 , 2 14 , 2 15 , and 2 16 scaled with the value 1.20 for the fractal dimension cfa which gives the best collapse of the plots 
in (a). 
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FIG. 4. 2D (a) Distributions P*(Ns\L) of the number of 3-blocks of mass N3 in a backbone of size L versus N3 for (from 
bottom to top) L = 16, 32, 64, 128, 256 and 512. The distributions exhibit a power-law regime with slope —2.25 =b 0.05 (b) 
Distributions scaled with the value 1.15 for the fractal dimension d% which gives the best collapse of the plots in (a). 
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x=N 3 /L da 

FIG. 5. 2D (a) Distributions P(Ns\L) of the number of top level 3-blocks of mass N3 in a backbone of size L versus 7V3 for 
(from top to bottom) L = 8, 16, 32, 64, and 128. The distributions exhibit a power- law regime with slope — 1.6 =b 0.1. (b) 
Distributions scaled with the value 1.15 for the fractal dimension ds which gives the best collapse of the plots in (a). 
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FIG. 6. 2D Distributions P(N 3 \L) of top level 3-blocks (filled symbols) and P*(N 3 \L) of all- level 3-blocks (unfilled symbols). 
While the slopes of the power law regimes of the two types of distributions are different, the finite-size-system cutoffs are 
essentially superimposed, consistent with the fractal dimension of the two types of distributions being equal. 
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FIG. 7. 2D (a) Distributions P* (Mz\L) of the number of 3-blocks of mass M3 in a backbone of size L versus M3 for from top 
to bottom) L — 16, 32, 64, 128, 256, and 512. In M3 we count not virtual bonds but all bonds in the 3-block. The distributions 
exhibit a power- law regime with slope —1.8 ±0.1 (b) Distributions scaled with the value 1.6 for the fractal dimension cfe which 
gives the best collapse of the plots in (a). 
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FIG. 8. 3D (a) Distributions P*(N3\N2) of the number of 3-blocks of mass N3 in a blob of size N2 versus N3 for (from 
bottom to top) N2 = 2 11 , 2 12 , and 2 13 . The distributions exhibit a power-law regime with slope —2.63 =b 0.05 (b) Distributions 
for N2 = 2 11 , 2 12 , and 2 13 scaled with the value 1.15 for the fractal dimension d% which gives the best collapse of the plots in 
(a). 
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x=N 3 /L° 3 



FIG. 9. 3D (a) Distributions P*(Ns\L) of the number of 3-blocks of mass N3 in a backbone of size L versus N3 for (from top 
to bottom) L = 8, 16, 32, 64, and 128. The distributions exhibit a power- law regime with slope —2.55 ±0.1. (b) Distributions 
scaled with the value 1.15 for the fractal dimension which gives the best collapse of the plots in (a). 
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x=N 3 /L da 

FIG. 10. 3D (a) Distributions P(Ns\L) of the number of top level 3-blocks of mass N3 in a backbone of size L versus N3 for 
(from top to bottom) L = 32, 64, 128, 256, and 512. The distributions exhibit a power- law regime with slope —2.0 ±0.1. (b) 
Distributions scaled with the value 1.15 for the fractal dimension ds which gives the best collapse of the plots in (a). 
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FIG. 11. 3D Distributions P(N 3 \L) of top level 3-blocks (filled symbols) and P*(N 3 \L) of all-level 3-blocks (unfilled symbols). 
While the slopes of the power law regimes of the two types of distributions are different, the finite-size- system cutoffs are 
essentially superimposed consistent, with the fractal dimension of the two types of distributions being equal. 
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FIG. 12. 3D (a) Distributions P*(Mz\L) of the number of 3-blocks of mass M3 in a backbone of size L versus M3 for (from 
top to bottom) L = 8, 16, 32, 64, and 128. In M3 we count not virtual bonds but all bonds in the 3-block. The distributions 
exhibit a power- law regime with slope — 1.87 =b 0.1. (b) Distributions scaled with the value 1.85 for the fractal dimension c/3 
which gives the best collapse of the plots in (a). 
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FIG. 13. (n*(N2)), the average number of 3-blocks in a blob of size N2 versus N2 for (a) 2D and (b) 3D. 
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